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1 .  INTRODUCTION 


Our  major  goal  In  the  work  reported  here  was  to  ascertain  whether 
the  electromagnetic  field  produced  by  the  AURORA  thick  target 
bremsstrahlung  could  be  modified  in  a  predictable  fashion  by  using  a 
parallel  plate  transmission  line.  If  this  prediction  can  be  shown,  then 
possibly  it  might  be  feasible  to  suggest  this  as  a  method  of  tailoring 
the  combined  electromagnetic  and  ionizing  radiation  so  that  it  more 
closely  resembles  the  characteristics  of  the  radiation  environment 
predicted  for  a  tactical  nuclear  weapon.  This  work  was  not  an  attempt 
to  design  a  tactical  electromagnetic  pulse  (EMP)  simulator,  but  rather 
was  directed  toward  gaining  insight  and  developing  tools  that  would  be 
of  value  in  designing  a  larger  test  facility,  one  that  could  be  used  to 
verify  EMP  source-region  coupling  calculations  and  possibly  even  test 
some  actual  tactical  systems.  The  experimental  measurements  reported  in 
this  paper  were  carried  out  in  the  fall  of  1977. 

2.  REQUIREMENT  FOR  MODIFICATION 

The  AURORA  Flash  X-Ray  Facility  at  the  Harry  Diamond  Laboratories 
was  designed  basically  to  study  high  intensity  transient  (nuclear) 
radiation  effects  on  electronics  TREE).  It  produces  high  intensity 
ionizing  radiation  over  a  relatively  small  volume.  The  resulting 
isodose  contours  produced  by  the  13-MeV-thick  target  bremsstrahlung  in 
the  metallically  shielded  AURORA  test  cell  are  shown  in  figure  1. 


The  components  necessary  to  simulate  tactical  nuclear  EMP  are  a  time 
varying  electric  (E-)  field,  a  time  varying  magnetic  (H-)  field,  a  time 
varying  Compton  electric  current,  and  time  varying  air  conductivity. 
Therefore,  at  first  look,  the  AURORA  test  cell,  which  incorporates  all 
these  components,  might  seem  to  be  an  ideal  candidate  for  the  simulation 
of  a  tactical  nuclear  EMP.  But  comparisons  of  measured  field  data*  and 
the  accepted  Defense  Nuclear  Agency  threat  criteria  indicate  serious 
differences  in  temporal  behavior  and  direction  and  magnitude  of  the 
various  components. 

We  can  gain  more  insight  into  the  problems  that  arise  by  considering 
the  possibility  of  using  the  unmodified  AURORA  as  a  tactical  nuclear  EMP 
simulator.  Let  us  simplify  matters  by  considering  the  wave  equations 
for  a  spatially  and  temporally  homogeneous  conducting  isotropic 
medium.  The  situation  inside  the  AURORA  test  cell  is  actually  not 
homogeneous,  and  the  medium,  the  ionized  air,  is  time  varying.  However, 
the  consideration  of  all  the  variables  would  require  the  use  of  such  a 
complicated  computer  code  that  essential  physics  would  tend  to  be 
obscured,  and  insight  into  the  actual  simulation  problem  would  be  lost. 

*Electric  field  data  discussed  later  in  this  section  and  unpublished 
magnetic  field  data  were  taken  by  Denis  Whittaker  of  the  AURORA  staff . 
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Figure  1.  AURORA  isodose  contours  and  transmission  line  in 
test  cell. 


Let  us  consider  the  wave  equations  inside  the  AURORA  test  cell. 


For  the  E-field,  we  have 


o2  l  3E  32e\  3Jc  Vp 

7  E  -  at  +  e°  3^j=  w +  r0  '  (1) 

where  o  is  the  air  conductivity,  Jc  is  the  Compton  current,  and  p  is  the 
charge  density  in  the  test  cell  ionized  air.  For  the  H-field,  we  have 


'2“  -  <•>  (•{?♦■.*)- 


=  -7  x  j 


If  we  had  no  Compton  current  and  resulting  air  conductivity,  the 
situation  would  be  greatly  simplified  because  the  above  two  complicated 


V 


wave  equations  would  become  the  two  relatively  simple  free  space  wave 
equations 


.  3  2p 

V2e  -  UoE0  —  *  0 
3t2 

and 

72h  -  u0e0  — -  =  0 
3t2 


(3a) 


(3b) 


and  our  simulation  problem  would  be  similar  to  the  exoatmospheric  EMP 
simulation. 

Now  consider  another  extreme:  suppose  the  conducting  current  were 
much  greater  than  the  displacement  current.  Then 


3e 

oe  »  e0 


(4) 


or  in  the  frequency  domain 


OE(<i>)  »  eoo)E(oj) 


(5) 


The  original 

wave  equations 

( 1 )  and 

( 2 )  now  become 

diffusion 

equations. 

_  2  3E 

72E  -  M0° 

3  J 

=  PO  -jf 

e0 

(6) 

and 


72H  -  p0a  ff  =  ~7  x  J 


(7) 
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where 


V2Ttfw0O 


(11) 


and 


E 

z 


Cle“Tx  +  C2*TX 


The  AURORA  test  cell  wall  is  a  ''perfect"  conductor;  therefore,  Ez  =  0 
alonq  the  test  cell  wall,  x  =  0,  and  we  have 


0  =  Ci  +  C2 


or 


Ci  =  -C2  =  (constant) 


Finally, 


Ez  =  (constant)  (e“x/5e“jx/6  -  e+x/6e+jx/6^  .  (12) 

Equation  (12),  which  applies  when  0E  >>  eQ (dE/dt) ,  is  an  expression 
that  can  be  employed  to  estimate  the  region  (distance),  6,  away  from  the 
conducting  AURORA  test  cell  walls  in  which  the  skin  effect  (in  the 
conducting  air)  short-circuits  the  E-field  parallel  to  the  AURORA  test 
cell  walls. 

Assume  that  o  =  10~3  mho/m;  then 


V2Tf4iM0"710"3 
1.125  (l04) 

VT 

For  f  -  108  Hz,  5  =  1.1  m;  for  f  =  109  Hz,  S  =  0.36  m. 


(13) 
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A  rough  estimate  of  the  horizontal  E-field  in  the  center  of  the 
AURORA  test  cell  when  oe  >>  eQ(3E/3t)  can  be  obtained  by  assuming  that 
the  E-field  builds  up  until  conduction  current  cancels  Compton  current? 
that  is, 

OE  =  -J  (14) 

c 

Actually,  there  would  probably  be  some  vertical  field  near  the  AURORA 
test  cell  floor,  ceiling,  and  walls  because  the  Compton  current  from  the 
AURORA  hot  spot  does  not  just  go  to  the  back  of  the  AURORA  test  cell, 
but  spreads  out,  and  a  portion  of  the  Compton  current  returns  to  the  hot 
spot  via  the  side  walls,  the  ceiling,  and  the  floor.  This  return  is 
illustrated  in  figure  2. 


Figure  2.  Sketch  of  average  Compton  current  trajectories. 


3.  MEASUREMENT  OF  AURORA  VERTICAL  ELECTRIC  FIELD 

When  the  AURORA  fires,  some  of  the  Compton  electrons  produced  by  the 
thick  target  bremsstrahlung  impinge  on  the  floor  and  generate  a  vertical 
E-field.  To  orient  ourselves,  we  need  to  know  the  magnitude  of  the 
vertical  E-field  on  the  test  cell  floor  even  though  it  may  not  be  the 
major  component  of  the  field  produced  by  the  AURORA  Compton  electron 
drivers . 


The  values  of  the  vertical  E-field  and  vertical  dE/dt  along  the 
AURORA  test  cell  floor,  5  m  from  the  front  test  cell  wall,  were  measured 
with  an  EGSG  PMD-1  E-field  sensor  and  with  a  parallel  grid  sensor 


suggested  by  Victor  van  Lint  (Mission  Research  Corp. ) .  The  grid  sensor 
was  connected  to  an  emitter  follower  with  a  high  impedance  input  so  that 
the  air  conductivity  did  not  influence  the  accuracy  of  the  E-field 

measurement.  Even  though  the  EGSG  sensor  was  designed  to  be  used  in  a 
vacuum,  its  values  agreed  to  within  20  percent  of  those  obtained  with 
the  van  Lint  sensor.  The  measurements  were  restricted  to  the  test  cell 
floor  because,  if  the  sensor  were  positioned  off  the  ground  plane,  the 
cables  to  the  sensors  would  perturb  the  field  that  the  sensor  was 
supposed  to  measure.  With  the  sensors  on  the  test  cell  floor,  the 

coaxial  cables  and  the  emitter  follower  could  be  buried  under  the  test 
cell  so  that  they  did  not  perturb  the  fields  in  the  test  cell.  Figure  1 
shows  the  AURORA  isodose  contours  and  the  position  of  our  first 
generation  vertical  E-field  enhancing  transmission  line. 

A  typical  vertical  E-field  measurement  at  a  distance  of  5  m  from  the 
front  of  the  test  cell  is  shown  in  figure  3.  We  need  to  increase  this 
field  if  we  want  our  simulator  to  produce  a  vertical  field  greater  than 

7  kV.  This  field  is  the  purpose  of  the  transmission  line.  The 

magnitude  of  the  vertical  E-field  is  probably  a  maximum  on  the  AURORA 
test  cell  floor.  If  the  vertical  E-field  were  measured  off  the  floor, 
intuition  and  symmetry  arguments  tell  us  that  the  field  would  decrease 
in  magnitude  and  then  change  direction  as  the  ceiling  were  approached. 
Van  Lint  is  designing  an  E-field  sensor  that  should  be  able  to  measure 
the  E-field  at  all  points  in  the  AURORA  test  cell. 


Figure  3.  Typical  vertical  electric-field  measurements  with 
Air  Force  Weapons  Laboraotry  sensor:  (A)  AURORA 
shot  2352,  (B)  fit  with  equation  (25),  and 


4.  TRANSMISSION  LINE  SIMULATOR 


The  use  of  a  parallel  plate  transmission  line  with  the  ionizing 
radiation  of  the  AURORA  to  simulate  the  electromagnetic  environment  that 
is  associated  with  a  tactical  nuclear  device  appears  to  be  a  basic 
approach  that  is  also  natural.  Consider  the  transmission  line  equations 
and  Maxwell's  equations: 


Transmission  line 
equations 


3v 

3x 


3i 

3x 


-  ("  ♦  *■  k)  1  ■ 

-  (°  ♦ c  h) v 


Field  equations 


V  x  E  =  -  M0  |f  , 


7  x  h  =  (a  +  e0 


Rut  R,  the  resistance  of  the  transmission  line  plate,  is  insignificant: 


3v 

3x 


f 


lj--(G  +  Ck)v  * 


Therefore,  for  the  transmission  line. 


32v  3v  32v 

3^7=  LC^  +  I^j^ 


and  for  Maxwell's  equations, 


3e  3  2E 

-  y0E0  +  JigEg 


at 


which  is  equation  (1)  without  the  driver  term.  If  H  ■  0,  the 
transmission  line  wave  equation  is  analogous  to  the  equation  derived 
from  Maxwell's  equations. 
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5.  LUMPED  PARAMETER  MODEL  OF  TRANSMISSION  LINE 

The  model  of  the  transmission  line  used  in  the  design  and  in  the 
calculated  predictions  for  our  experiments  is  a  lumped  parameter 
transmission  line  consisting  of  12  sections.  Figure  4  shows  the  basic 
section  that  describes  the  situation  in  which  the  transmission  line  is 
pulsed,  but  the  AURORA  is  not. 


Figure  4.  Pulsed  transmission  line. 


The  value  of  the  capacitance  in  farads  is  determined  by  the 
expression 


c  =  e0A/d 


where  is  the  permittivity  constant,  A  is  the  area  of  the  section,  and 
d  is  the  separation  of  the  transmission  line. 

The  value  of  the  inductance,  L,  was  determined  by  employing  equation 
(15)  for  the  impedance  of  the  transmission  line  in  terms  of  L  and  C: 


The  impedance  of  the  transmission  line  was  determined  by  the  technique 
of  time  domain  ref lectometry.  That  is  to  say,  the  reflection  of  the 
transmission  line  was  found  to  be  insignificant  when  it  was  terminated 
with  a  90-fi  resistor.  We  could  then  use  equation  (15)  to  obtain  the 
value  of  the  inductance,  or  L  *  (90)^C.  With  L  and  C  determined,  we 
were  in  the  position  to  introduce  the  effect  of  time  varying  air 
conductivity  into  a  lumped  parameter  section.  The  transmission  line 
section  was  modified  to  incorporate  the  effect  of  a  time  varying 
resistor  in  parallel  with  the  capacitor  as  shown  in  figure  5. 
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L/2  L/2 


Figure  5.  Modified  transmission  line. 


The  time  varying  resistor,  R^,  is  assumed  to  vary  as 
sin-2  {itt/[t  +  2  exp  (-In  2  T01)]}  rad  (Si)/s 


where 

T  =  t/t* 


t*  =  time  of  peak, 
a  =  1.77  t*/doseJ  1  • 9  , 

A*  *  peak  dose  rate  . 


This  curve  seems  to  fit  a  variety  of  AURORA  pulse  shapes  quite  well.  It 
is  a  smooth  curve  and  lacks  the  bumps  and  occasional  early  time  foot 
seen  in  some  shots.  Digitized  shot  data  would  be  better,  but  this  fit 
seems  adequate.  In  a  few  of  the  calculations  not  compared  with  measured 
data,  the  simple  dependence  cos-2  [ir(T  -  TQ/2)]/T0  (with  TQ  =  3  x  io“7 
s)  was  used. 

Air  chemistry  and  air  conductivity  calculations  can  be 
complicated.  However,  the  relatively  slow  rise  time  of  the  AURORA  pulse 
(10  shakes)  and  the  relatively  fast  attachement  rate  of  electrons  to  the 
oxygen  component  of  air  at  sea  level  ( 1  shake)  allow  us  to  obtain  a 
simple  relationship  between  the  AURORA  gamma  flux  and  the  resulting  air 
conductivity. 

In  brief,  it  can  be  shown  that,  in  sea  level  air,  the  ionization 
rate  of  a  directed  gamma  flux  is  roughly  y  =  2  x  109FY,  where  y  is  the 
ionization  rate  in  ion  pairs/cm3s,  and  is  the  directed  gamma  flux  in 
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roentgens/s.  The  attachment  rate  of  electrons  in  sea  level  air 
approximately  equals  10~®  s”'-.  Since  the  AURORA  gamma  flux  varies 
slowly  compared  with  the  electron  attachment  rate,  the  electron  density, 
Ne,  can  be  estimated  by 


electronsX  I  .  20F  .  (16) 

cm3  /  “  Y 

In  a  collision  dominated  plasma  such  as  partially  ionized  air  at  sea 
level,  the  air  conductivity,  a,  is  given  by 


a 


e2N 


mv 


eN  y 
e 


(17) 


where  m  is  the  electron  mass,  v  is  the  collision  frequency  between  free 
electrons  and  air  molecules,  and  y  is  the  electron  mobility.  Equations 
(16)  and  (17)  can  be  combined  to  yield* 


ae 


5  x  10-14N 

e 


(18) 


=  0.5  x  10~12F  . 

Y 

The  gamma  flux  can  be  estimated  from  the  total  gamma  dose  and  the  AURORA 
pulse  shape.  The  gamma  dose  is  measured  with  thermoluminescent 
dosimeters,  and  the  directed  flux  or  the  dose  rate  is  then  calculated  by 
using  the  measured  AURORA  pulse  shape.  At  a  later  time,  after  a  number 
of  our  experiments  involving  the  effect  of  humidity  on  the  electron 
mobility  have  been  digested  and  analyzed,  we  will  allow  the  electron 
mobility  to  be  a  function  of  the  E-field  and  the  air  humidity.  In  our 
present  analysis,  we  use  the  simple  relation  given  in  equation  (18).  If 
the  conductivity,  o^(T),  of  a  transmission  line  section  is  known,  we  can 
obtain  the  resistance,  R^(T).  Let  the  area  of  the  transmission  line 
section  be  A  and  the  distance  between  the  upper  and  lower  conductors  be 
d ;  then 


Ri(T) 


d 

0.(T)A  ' 


*The  value  0.5  x  20“12  F  is  obtained  from  Dean  May's  analysis  of  our 
" pie  pan "  air  conductivity  measurements.  Dean  L.  May,  Electron  Mobility 
and  Electron  Attachment  Rate  in  an  Electromagnetic  Pulsed  Air 
Environment,  HDL  Student  Technical  Symposium  (17-18  August  1977),  87. 
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and  finally 


-_y  | 

R.  (T)  =  - - —  sin-2  {irr/[T  +  2  exp  (-In  2  T  )]}  ,  (19) 

1  O  .A  I 

maxi  I 

where  R^(T)  is  the  time  varying  resistance  across  transmission  line 

section  i.  Figure  5  shows  the  transmission  line  section  with  the  time 
varying  resistor  across  the  capacitor  to  approximate  the  effect  of  the 
time  varying  conductivity. 

Another  physical  effect  has  to  be  incorporated  in  our  model:  the 
effect  of  the  AURORA  generated  vertical  E-field  must  be  superimposed  on 
our  transmission  line  section.  It  can  be  superimposed  by  inserting  a 
current  source  in  parallel  with  the  capacitor  and  the  time  varying 
resistor.  The  modification  is  shown  in  figure  6.  As  seen  in  figure  6,  *i 

the  current  source  consists  of  a  conduction  current  and  a  displacement 
current.  I 

L/2  L/2 


We  previously  described  our  measurement  of  the  values  of  the 
vertical  E-field  and  dE/dt  produced  in  the  "empty”  AURORA  test  cell  (the 
AURORA  test  cell  without  the  transmission  line)  when  the  AURORA  is 
fired.  The  E- fie Id  used  in  our  current  source  is  shown  in  figure  2. 

We  can  now  write  the  series  of  voltage  and  current  first-order 
differential  equations  that  describe  the  transmission  line,  but  first  I 

let  us  consider  how  the  first  transmission  line  section  is  connected  to 
the  pulser:  A  low  inductance  2-MF  capacitor  is  discharged  into  the 
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mummed 


transmission  line  as  shown  in  figure  7.  To  be  consistent  with  our 
Runge-Kutta  subroutine,  we  label  the  voltages  with  odd  numbers  and  the 
currents  with  even  numbers*  In  other  words,  the  voltage  at  node  N  is 
designated  V(2N  -  1),  and  the  current  flowing  in  the  line  into  node  N  is 
labeled  I(2N).  The  advantages  of  this  method  of  notation  are  clear  when 
we  consider  the  actual  node  equations  to  be  solved  by  the  Runge-Kutta 
method . 


V(1)  RO 


M  +  eQ  dE_ 

R1  +  dT  A 


TRIGGER 


HIGH 

VOLTAGE 

POWER 

SUPPLY 


— WV-J” 


Figure  7.  Discharge  into  transmission  line. 

Let  us  consider  the  partial  circuit  shown  in  figure 
first  current  and  voltage  equations  are 

7.  Three 

of  the 

CO  =  -1(2)  , 

dT 

(20) 

LI  —  ■j|)-  =  V(1,  -  V(  3 )  -  RO  x  1(2) 

t 

(21) 

dV { 3 )  V(3)  T...  ,  ABC 

C1  dT  -  "  R1(T)  +  I(2>  '  I(4)  +  R1(T)  + 

( DABC )  , 

(22) 

where 

The  values  of  E  and  dE/dT  in  equations  (23)  and  (24)  are  given  by 
equations  (25)  (multiexponential  field)  and  (26),  respectively.  The 
values  of  the  constants  in  equations  (25)  and  (26)  were  obtained  by  the 
curve  fit  shown  in  figure  3.  In  our  first  application  of  the  code,  we 
assume  that  E  and  dE/dT  are  constant  along  the  transmission  line.  we 
do,  however,  allow  the  value  of  crfT)  or  R^(T)  to  be  determined  by  the 
radiation  dose  along  the  transmission  line. 


E(T) 


( A1  -  A2 ) ( 1  +  nT)e“nT  +  (A2  +  A3)(1  +  6T)e 


-6t 


-aT  .  -VT 

-  A1(1  +  aT)e  -  A3 ( 1  +  YT)e 


—6  '  f T— T* 1 

+  A4  U(T  -  T  '  )  {  [  1  +  6'(T  -  T’)]e 


(25) 


-a • ( T-T ' ) i 

-  [1  +  a'(T  -  T*  )]e  '  '} 


=  A 1  a2Te-aT  +  A3  Y2Te-YT 
dT 


-  ( Al  -  A2)n2Te*nT  -  (A2  +  A3)62Te“0T 


(  „  o  - & ' ( T-T  * )  ( 26) 

-  A4  U(T  -  T')(.(3')2(T  -  T')e 

-  (a* )2(T  -  T* )e-a,(T’T’) j 


Figure  8  shows  the  circuit  for  the  termination  end  of  the  transmission 
line. 
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Figure  8.  Termination  of  transmission  line. 


The  current  sources  I.  =  E-d/R^  +  e^ACdE/dT)  vary  as  a  function  of 
position  along  the  transmission  line  because  R^  varies  as  a  function  of 
radiation  dose.  In  a  more  accurate  calculation,  E  and  dE/dT  would  most 
likely  vary.  Also,  R^  would  be  a  function  of  dose  and  E-field  because 
mobility  varies  as  a  function  of  E.  The  effect  of  the  variation  of  A  in 
the  tapered  sections  of  the  transmission  line  would  affect  also  the 
value  of  the  displacement  current,  £QA(dE/dT). 

The  24  differential  equations  that  describe  the  lumped  parameter 
model  of  the  transmission  line  are  given  in  table  I  (p.  22).  The  table 
is  actually  a  portion  of  the  Runge-Kutta  FORTRAN  subroutine. 


6.  TRANSMISSION  LINE  EXPERIMENTS 

6.1  Physical  Description  of  Transmission  Line 

The  transmission  line  was  fired  with  three  different  pulser 
configurations  (fig.  9  to  11). 

Figure  12  shows  a  schematic  view  of  the  transmission  line,  and 
an  alteration  of  the  transmission  line  configuration  was  to  shield  the 
tapered  ends  of  the  transmission  line  with  a  lead  brick  wall.  The 
purpose  of  the  lead  shield  was  to  cut  down  the  ionization  of  the  air  in 
the  tapered  sections  of  the  line  and  to  deliver  as  much  E-field  and 
energy  to  the  nontapered  sections  of  the  transmission  line  as 
possible.  Various  configurations  of  the  transmission  line  with  and 
without  lead  shields  are  shown  in  figure  13. 
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TABLE  I.  EQUATIONS  TO  SOLVE  TRANSMISSION  LINE  TERMINATED  IN 
90-ft  RESISTANCE 


F  C 1 )  =-Y(2)/C0 

F(2)  =  Y ( ll/Ll-Y (3J/L1  - (R  0/Ll ) *Y  12  ) 

F ( 3 )  =  -Y(3)/(  R 1  *C  1  )  +Y(2)/C1  -Y(4)/C1  *  A6C/(R1*C1)  4  DABC 

F  ( <> )  =  Y(3)/L2  -  Y  {  5  )  /  L2 

F(5)  *  -Y (5 )7  ( C2*R2 )  +  Y(4)/C2  -  YC61/C2  +  ABC/CC2-R2')  +  DABC 

F (6)  =  Y(5)/L3  -  Y17J/L3 

F<7)  =  -  Y(7)/(C3-R3)  +  Y16J/C3  -  Y(8)/C3  ♦  A3C/(C3*R3)  ♦  DABC 

F{'61  =  Y  i  7) /  L4  -  Y  (9)7 LA 

F ( 9 )  =  -Y(9)/(C4=R4)  +  Y{8)/C4  -  Y(10)/C4  +  ABC/IC4*R4)  4  DABC 

F ( 1 0 1  =  Y { 9 ) / L5  -  Y  C 1 2  )/L5 

F  (11  )  =  -Y(ilT/(C5«R5)'  +  Y  { 10J/C5  -Y(12)/C5  4  A3C/IC5*R5)4  DABC 
F  1 1 2  1  =  Y  ( 1 1  J  /L  6  -  YU21/L6 

F  (  1 3  )  =  -Y  { 1  3  ) /(C6-R6  )  +  Y1121/C6  -VUM/C6  4  ABC/tC6*R6)  4  DABC 
F ( 1 A )  =  Y ( 13 ) / L7  -  Y(15)/L7 

F ( 1 5 )  =  -Y ( 1 5 ) / ( C7*R7 )  4  Y1141/C7  -  Y1161/C7  4  ABC/(C7*R7) 

1 _ 4 _ DABC  _ _ _ 

F  {  1 6  )  =  YU5J/L8  -  YTl7)/L8 

F  C 17 1  =  -Y ( 1 7 ) / { C8  ®R  8 1  4  Y(16)/C8  -  YC18)/C8  .+  ABC/CC8*R8) 

1_  4  JABC _ 

">(181  =  YC171/L9  -YC191/L9 

F ( 19 1  =  -Y 119 J/IC9-R9)  +  Y118J/C9  -Y(20)/C9  4  ABC/CC9-R9) 

1  4  DABC  _  _ ___  _  _  _ 

F ( 20 1  =  Y(19}/L10-Y(21  )/L10 

F ( 21 )  *  -Y (21 )/«ClC«R10)4Y l 20 ) /C 10-Y l 22) /C 10  4  ABC/ ( R1 0*C 10  ) 

1  4  DABC  _  _ 

F ( 22 1  =  Y ( 21 )  /L i 1 -Y I  23  )  /U  1 

F  ( 23  )  =  -YI23)/(C11*RM  )  >Y ( 22  )/C 1 1-Y ( 24  ) /C 1 1  ♦  ABC/ (  C 1 1~R  1 1 ) 

1  4  DA3C  _  _ 

'  F12M  =  Y  (  23)  /LI  2  - Y  (  2 A  )  «  l  R 1 2 7 L 1 2  ) 


Note:  Let  N  refer  to  node  N, 
dV(N)/dT  =  F (2N-1) , 

V(N)  =  Y (2N-1) 
dl (N)/dT  =  F(2N) , 

I(N)  =  Y  (2N)  , 

where  I(N)  is  the  current  flowing  into  node  N. 
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Figure  9.  Pulser  configuration  I. 


Figure  11.  Pulser  configuration  III. 
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6.2  Experimental  Measuring  Techniques 

A  number  of  parameters  were  measured  during  each  AURORA  shot 
involving  the  transmission  line.  Two  E-field  sensors  were  placed  in  the 
parallel  section  of  the  transmision  line.  One  sensor  was  the  Air  Force 
Weapons  Laboratory  (,\FWL)  EG&G  PMD-1  sensor  designed  by  Carl  Baum;  the 
other  was  a  parallel  plate  sensor  that  fed  into  a  high  impedance  input 
emitter  follower  to  increase  the  sensor  time  constant.  Its  general 
design  was  suggested  by  van  Lint.  At  the  center  of  the  transmission 
line,  a  Moebius  loop  was  installed  to  measure  the  H-fleld.  Current  into 
the  transmission  line  and  current  at  the  matched  termination  were 
measured  with  ionizing  radiation  insensitive  Adelco  current  probes. 


7.  COMPARISON  OF  EXPERIMENT  AND  CALCULATION 

Four  electromagnetic  parameters  were  measured  with  each  simultaneous 
firing  of  the  AURORA  for  the  transmission  line.  These  experimental 
parameters  were  the  E-field  at  the  center  of  the  transmission  line,  the 
H-field  at  the  '•inter  of  the  transmission  line,  the  current  at  the 
termination  of  the  transmission  line,  and  the  current  at  the  input  of 
the  transmission  line.  We  measured  also  the  total  ionizing  radiation 
dose  at  various  points  along  the  transmission  line.  For  this  report,  we 
only  digitized  our  data  and  corrected  them  with  a  measured  transfer 
function.  We  did  not  apply  all  possible  corrections  and  refinements  to 
our  theory.  Air  conductivity  refinements  such  as  humidity  and  radiation 
dose  variations  from  shot  to  shot  were  left  for  later  consideration. 
Also,  a  few  systematic  errors  needed  to  be  resolved. 


7. 1  Measured  Electric  Field  and  Calculated  Transmission  Line 
Voltage 

7.1.1  Configuration  I 

The  E-field  in  the  transmission  line  was  measured  with  the 
AFWL  EG&G  E-field  sensor  and  with  an  E-field  sensor  designed  by  van 
Lint.  Since  the  parallel  transmssion  line  separation  was  0.4  m,  the 
transmission  line  voltage  is  given  by  V  =  0.4E.  Figure  14  compares  the 
calculated  and  measured  E-fields  across  the  untapered  portion  of  the 
transmission  line  for  impedance  matching  configuration  I  when  the 
transmission  line  pulser  was  not  pulsed;  the  voltage  across  the 
transmission  line  section  was  caused  only  by  the  AURORA  generated 
electromagnetic  field  in  the  AURORA  test  cell. 

Figure  15  compares  the  calculated  and  measured  E-fields  for 
experimental  configuration  I  when  the  pulser  capacitor  was  charged  to 
6000  V  and  discharged  a  few  micrseconds  before  the  AURORA  was  fired. 
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TIME  1*10  8  * 

Figure  14.  Electric  field,  transmission  line  pulser  not 
pulsed,  configuration  I. 


Figure  15.  Electric  field,  transmission  line  pulser  pulsed, 
configuration  I. 


7.1.2  Configuration  II 


Figure  16  shows  the  calculated  and  measured  E-field  across 
the  untapered  portion  of  the  transmission  line  for  impedance  matching 
configuration  II.  The  transmission  line  pulser  was  not  pulsed  so  that 
the  voltage  appearing  across  the  transmission  line  was  caused  only  by 
the  AURORA  generated  electromagnetic  field  in  the  AURORA  test  cell. 


TIME  l»108  .1 

Figure  16.  Electric  field,  transmission  line  pulser  not 
pulsed,  configuration  II. 


Figure  17  compares  the  calculated  and  measured  E-fields  for 
experimental  configuration  II,  except  that  the  transmission  line 
capacitor  was  charged  to  4000  V  and  discharged  a  few  microseconds  before 
the  AURORA  was  fired. 

Figure  18  shows  a  purely  calculated  result.  The  transmission 
line  was  assumed  to  be  in  impedance  matching  configuration  II,  and  the 
electromagnetic  field  in  the  AURORA  test  cell  was  assumed  to  be  zero, 
but  the  ionizing  radiation  was  assumed  to  be  present.  The  transmission 
line  capacitor  was  assumed  to  be  charged  to  4000  V,  and  the  pulser  was 
assumed  to  be  discharged  about  100  ns  before  the  AURORA  was  fired. 


7.1.3  Conf iguration  III 


When  the  transmission  line  tapered  sections  were  not  shielded 
with  lead,  the  unreliability  of  both  the  trigger  circuit  and  the 
accompanying  trigatron  interfered  with  measurements  for  impedance 
matching  configuration  III.  To  not  exceed  the  dynamic  range  of  the  E- 
field  measuring  electronics,  the  transmission  line  capacitor  was  not 
charged  beyond  2000  V.  Unfortunately,  the  triggering  circuit  and  the 
trigatron  could  not  reliably  operate  at  such  a  low  voltage.  However, 
two  interesting  calculational  results  correspond  to  configuration  III. 

Figure  19  shows  the  calculated  behavior  of  the  transmission 
line  when  the  measured  E-field  in  the  AURORA  test  cell  is  taken  into 
consideration.  The  large  late  time  spike  occurs  whenever  the  impedance 
matching  network  does  not  limit  the  current  that  can  flow  into  the 
transmission  line.  A  possible  explanation  for  the  spike  is  that  the 
ionizing  radiation  short-circuits  the  transmission  line,  and  the 
resulting  large  current  produces  a  large  H-field  and  an  accompanying 
large  magnetic  energy  in  the  line.  When  the  ionizing  radiation 
disappears,  some  of  the  magnetic  energy  is  converted  into  electrical 
energy,  and  a  large  voltage  is  produced  across  the  transmission  line. 

The  results  shown  in  figure  20  are  similar  to  the  results 
shown  in  figure  19,  except  that  the  E-field  produced  by  the  AURORA  in 
the  test  cell  is  not  included  in  the  figure  20  calculation. 


Figure  19.  Calculated  electric  field  of  transmission  line, 
configuration  III. 
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Figure  20.  Calculated  electric  field,  transmission  line, 
configuration  III  (AURORA  electric  field 
not  included) . 


7.2  Experimental  Measurements  with  Lead  Shielding  in  Place 

Figure  21  compares  the  calculated  and  measured  E-fields  for 
impedance  matching  configuration  I;  the  transmission  line  was  pulsed 
with  6000  v,  and  the  tapered  section  of  the  transmission  line  near  the 
pulser  was  shielded  with  lead.  Figure  22  compares  the  calculated  and 
measured  E-fields  for  a  configuration  II  impedance  match;  the 
transmission  line  was  pulsed  4000  V,  and  the  tapered  section  near  the 
pulser  was  shielded  with  lead.  Figure  17  shows  the  equivalent 
unshielded  result.  Figure  23  compares  the  calculated  and  measured  E- 
fields  for  a  configuration  II  impedance  match.  The  AURORA  and  the 
transmission  line  pulser  were  fired  at  approximately  the  same  time.  For 
this  measurement,  the  tapered  transmission  line  transition  sections  both 
near  the  pulser  and  next  to  the  termination  were  shielded  with  2- in. 
(5.08-cm)  lead  shield  walls. 
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Figure  21.  Calculated  and  measured  electric  fields,  transmission 
line,  configuration  I. 


Figure  22.  Calculated  and  measured  electric  fields, 
transmission  line,  configuration  II. 


TIME  I*  10  *  l) 

Figure  23.  Calculated  and  measured  electric  fields, 
transmission  line,  configuration  II 
(AURORA  and  pulser  fired  at  approximately 
same  time) . 

7.3  Measured  and  Calculated  Input  Currents 

The  lumped  parameter  code  also  can  be  used  to  calculate  the 
current  into  the  transmission  line  at  the  points  indicated  in  figures  10 
and  11. 

Figure  24  compares  the  calculated  and  measured  input  currents 
to  the  transmission  line  for  input  impedance  matching  configuration 
II.  Figure  25  compares  the  input  current  and  the  calculated  current  for 
input  configuration  III.  Figure  24  corresponds  to  the  unshielded  shot, 
and  figure  25  corresponds  to  the  shot  in  which  both  tapered  ends  of  the 
transmission  line  are  shielded  with  lead.  The  model  predicts  the  input 
current  to  within  50  percent. 

7.4  Measured  and  Calculated  Termination  Currents 

Figure  26  compares  the  calculated  and  measured  termination 
currents  for  transmission  line  pulser  impedance  matching  configuration  I 
with  the  transmission  line  pulser  capacitor  charged  to  6  kV.  Figure  27 
compares  the  calculated  and  measured  termination  currents  for  impedance 
matching  termination  configuration  I,  but  with  the  transmission  line 
pulser  capacitor  charged  to  0  V. 
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Figure  24.  Calculated  and  measured  input  currents 


transmission  line,  configuration  II. 


Figure  25.  Calculated  and  measured  input  currents 
transmission  line,  configuration  III. 


Figure  27.  Calculated  and  measured  termination  currents, 

transmission  line,  configuration  X  (capacitor  0  V) . 


Figure  28  shows  the  calculated  and  measured  termination 
currents  for  pulser  impedance  matching  configuration  II  with  the 
transmission  line  pulser  capacitor  charged  to  4  kV.  Figure  29  shows  the 
calculated  and  measured  termination  currents  for  impedance  matching 
configuration  II  with  the  capacitor  charged  to  0  V. 

Figure  30  compares  the  calculated  and  measured  termination 
currents  for  pulser  impedance  matching  configuration  III.  Figure  30 
corresponds  to  an  experimental  setup  with  the  tapered  section  of  the 
transmission  line  near  the  pulser  shielded  with  a  5.08-cm  lead  wall. 
The  reason  that  we  were  able  to  obtain  these  results  is  that  the  pulser 
was  charged  to  3  kV.  This  high  voltage  saturated  the  E-field  sensor, 
but  allowed  us  to  synchronize  the  transmission  line  pulse  with  the 
AURORA  pulse.  The  pattern  emerging  from  figures  26  to  30  is  that  our 
lumped  parameter  model  of  the  transmission  line  is  a  qualitative  and 
even  semiquantitative  prediction  of  the  actual  measured  termination 
current. 


Figure  28,  Calculated  and  measured  termination  currents, 

transmission  line,  configuration  II  (capacitor  4  kV) . 
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7.5  Magnetic  Field  Measurement 

The  Moebius  loop  H-field  sensor  was  oriented  so  that  it  would 
measure  the  H-field  produced  by  the  current  flowing  along  the 
transmission  line.  Unfortunately,  only  one  Moebius  loop  was  positioned 
in  the  transmission  line,  and  it  was  at  one  half  the  distance  between 
the  transmission  line  pulser  and  the  transmission  line  termination. 
Although  it  probably  would  have  been  more  informative  if  we  had  had  many 
H-field  sensors  distributed  along  the  transmission  line,  the  number  of 
data  channels  and  oscilloscopes  limited  their  placement. 

According  to  our  lumped  parameter  model ,  the  shape  of  the 
voltage  pulse  does  not  vary  appreciably  from  node  to  node;  however,  the 
shape  of  the  current  pulse  does  vary  appreciably  from  the  transmission 
line  pulser  to  the  termination.  The  explanation  may  be  found  in  terms 
of  the  general  transmission  line  equation: 


where  R,  L,  C,  and  G  are  shown  in  figure  31.  Usually,  R  and  G  are  very 
small,  and  the  transmission  line  equation  for  a  coaxial  line  then 
reduces  to  the  familiar  form 


32v  3v2 

—  *  1/2  — r  * 

3x2  3t2 

32i  T_  32v 

3x2  3t2 


In  our  experiment,  G  is  a  current  sink  that  drains  off  current  as  we 
progress  along  the  line.  We  have  not  only  a  current  sink,  but  also  a 
current  source  in  parallel  with  the  sink  (sect.  5). 
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Figure  31.  Schematic  of  transmission  line. 


Another  more  succinct  drawing  of  a  section  of  a  lumped 
parameter  transmission  line  is  shown  in  figure  32.  In  this  line. 


3 

Z  =  R  +  L  t—  (ohm/m) 


and 

Y  =  G  +  C  r—  (mho/m) 


Figure  32  depicts  an  approximate  model.  There  are  ways  of  obtaining  an 
exact  model  for  a  section  of  line  at  a  single  frequency,  but  we  are 
working  in  the  time  domain. 


ZAX  ZAX 


Figure  32.  Lumped  parameter  transmission  line. 
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A  model  with  the  above  approximate  sections  is  a  reasonably 
good  representation  of  a  transmission  line,  provided  that  the  voltage 
drop  between  points  1  and  2  is  small  compared  with  the  voltage  across 
the  line  at  either  end  of  the  section  and  that  the  loss  in  current 
through  the  shunt  admittance  is  small  compared  with  the  total  current 
flowing  in  the  line.  That  is. 


V2  -  Vi  <<  Vi 


and 


X2  -  Ii  «  Ii 


The  first  voltage  inequality  was  met  in  our  calculations  and 
measurements.  However,  the  current  inequality  was  met  only 

marginally.  Our  H-field  measurements  were  intrinsically  related  to  the 
current  in  the  line.  From  a  theoretical  point  of  view,  we  most  probably 
can  improve  our  model  by  increasing  the  number  of  sections  in  the  lumped 
parameter  model.  Ultimately,  as  the  number  of  sections  is  increased, 
both  inequalities  will  be  met.  Hiis  matter  can  be  investigated  by 
computer  experiments  in  which  the  number  of  sections  is  increased. 

Consider  the  transmission  line  in  the  case  of  a  pulser 

Impedance  match  of  type  I  with  the  pulser  capacitor  charged  to  6  kV. 
Figures  33  and  34  show  the  voltage  and  the  current  at  node  2,  and 
figures  35  and  36  show  the  voltage  and  the  current  at  termination  node 
12.  The  voltage  pulses  at  nodes  2  and  12  are  quite  similar,  but  the 
current  pulses  are  quite  different.  Figure  37  plots  the  current  maximum 
or  minimum  ■"nring  the  AURORA  pulse  as  a  function  of  node  number  or 

position  along  the  line.  Figure  37  emphasizes  the  variation  of  current 
pulse  shape  as  a  function  of  position  along  the  line. 

Figure  38  compares  the  variation  of  the  current  inferred  from 
the  Moebius  loop  and  the  calculated  current  flowing  into  node  7.  In 
general,  we  can  say  that  figure  38  shows  fairly  good  agreement  between 
calculation  and  measurement.  However,  a  more  definitive  comparison  with 
H-fields  along  the  transmission  line  and  our  lumped  parameter  theory 
would  require  a  number  of  Moebius  loops  spaced  along  the  transmission 
line  rather  than  just  one  loop  halfway  between  the  pulser  and  the 

|  termination. 

|  Figures  39  and  40  are  calculated  currents  for  nodes  6  and  7  for 

the  capacitor  charged  to  0  V  and  the  pulser  matching  circuit  of 

configuration  II.  The  current  changes  sign  near  the  center  of  the  line 
for  this  configuration  and  capacitor  voltage;  thus,  it  is  a  poor 
X  position  for  our  single  Moebius  loop. 
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8.  CONCLUSIONS 

This  report  describes  a  series  of  experiments  in  which  a  parallel 
plate  transmission  line  was  used  to  augment  the  electromagnetic  field  in 
the  AURORA  test  cell  during  an  AURORA  shot.  A  simple  lumped  parameter 
model  of  the  transmission  line  was  modified  with  time  varying  resistors 
and  with  parallel  current  sources.  This  model  is  a  good  tool  to  use  to 
understand  the  many  aspects  of  the  behavior  of  the  line  when  subjected 
to  the  time  varying  ionizing  radiation  produced  by  the  AURORA  firing.  A 
more  sophisticated  calculation  is  being  developed  for  the  time  varying 
resistances  across  the  transmission  line  during  an  AURORA  shot.  When 
results  obtained  from  that  experiment  are  compared  with  the  finally 
analyzed  measurements,  the  agreement  may  very  well  become  closer.  We 
are  confident  that  our  transmission  line  theory  gives  a  qualitative, 
even  semiquantitative  description  of  the  superposition  of  a  vertical  E- 
field  upon  the  conducting  ionized  air  by  the  transmission  line.  The 
next  task  is  to  use  this  simple  transmission  line  model  to  determine 
whether  or  not  to  modify  the  AURORA  test  cell  environment  so  that  it 
more  closely  resembles  that  associated  with  the  actual  tactical  source 
region. 
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